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Abstract 

The concept of cross diffusion is applied to some biological systems. 
The conditions for persistence and Turing instability in the presence 
of cross diffusion are derived. Many examples including: predator- 
prey, epidemics (with and without delay), hawk-dove-retaliate and 
prisoner's dilemma games are given. 

1 Introduction 

Cross diffusion is the diffusion of one type of species due to the presence of 
another [Okubo 1980, Chattopadhyay and Chatterjee 2001]. This phenomena 
is abundant in nature e.g. predator-prey systems where the predator diffuses 
towards the regions where the prey is more abundant. On the other hand, the 
prey tries to avoid predators by diffusing away from them. Another area of 



application is in epidemics where susceptible individuals try to avoid infected 
individuals. 

Our aim is to generalize some biological systems to the case of the presence 
of cross diffusion. We study the persistence and Turing instability of the 
generalized models. Several examples are given e.g. predator-prey system, 
epidemic model, hawk-dove-retaliate game and prisoner's dilemma game. 

The paper is organized as follows: In section 2 the conditions for Turing 
instability in the presence of cross diffusion are derived then applied to a 
predator-prey system. Necessary conditions for persistence of a predator 
prey system in the presence of cross diffusion are derived. In section 3, cross 
diffusion is introduced in an epidemic model. Its effect on wave solution 
with and without delay is derived. In section 4, the effect of cross diffusion 
on evolutionarily stable strategy is studied. Hawk-dove-retaliate game is 
given as an example. Cross diffusion is applied to the dynamics of learning 
in multiagent systems in section 5. Some conclusions are summarized in 
section 6. 



2 Cross diffusion in a predator-prey model 

Cross diffusion is expected to be relevant to predator-prey systems where 
the predator diffuses towards the regions where the prey is more abundant. 
On the other hand, the prey tries to avoid predators by diffusing away from 
them. Cross diffusion term expresses the population flux of one species due 
to the presence of another species. The value of the cross diffusion may be 
positive, negative or zero. 

Positive cross diffusion means that one type of species tends to move in 
the direction of lower density of the other type and vice versa. Hence the 
cross diffusion term is positive for the prey and negative for the predator. 
Generalizing Lotka-Volterra predator-prey model (with prey saturation) to 
include cross diffusion, one gets 

% = uMi - -u 2 )+D 1 ^ + £(D' 12 ( Ul )^), 

% = u 2 (-a 2 + Ul ) + D 2 ^ - £(D> 21 (u 2 ) 9 i±), { ) 

where D' 12 (u\) is the density dependent diffusion coefficient of cross diffusion 
such that D' 12 (ui) — > as u\ — > 0. Following Chattopadhyay and Chatterjee 
[Chattopadhyay and Chatterjee 2001] we assume that 

D' 12 = Di 2 ui/(e + ui} « D 12 (l - eju x + ...), 
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where e is very small hence 



D' 12 ks D 12 Vui >> e. 
Similarly for D 2 i thus we get 



^ = - - « 2 ) + + £>i 2 ^, r2 x 

^ = M2 (-a 2 + Ml ) + ^fr-^i&, 1 ] 

where 1*1(1x2) is the prey (predator) density, and all the constants cci, a 2 , Pi, 
Di, D 2 , D i2 , D 2 i are nonnegative. We have used —D 2 i since predators move 
to regions with higher prey density. The value a± is the net growth rate of 
the prey population in the absence of predators, while a 2 is the net death 
rate of predators due to the absence of prey. The term P1U1 corresponds to 
the competition within the prey species. 

The existence of solutions of systems similar to (2) has been studied in 
[Chen et al 2003 and Le 2003]. 

This system has a unique homogeneous coexistence solution 

u* = a 2 , u 2 = a.\ — a 2 p±, ol\ — a 2 f3\ > 0. (3) 

It is straightforward to see that the solution (3) is asymptotically stable in 
the case of no diffusion (D 1 = 0, D 2 — 0, _D 2 i = 0, D 12 = 0). 

In Turing instability [Okubo 1980], it is required to see whether diffusion 
terms can destabilize the steady state (2). We derive Turing instability in 
the presence of cross diffusion: Consider two interacting species with different 
diffusion coefficients 

— = f 1 ( Ul ,u 2 ) + D 1 — + D 12 — , (4) 
du 2 . , d 2 u 2 d 2 u\ 



= f 2 ( Ul ,u 2 ) + D 2 —^ + D 



at ox 1 ox 1 

Assume that {u\,u^} are the spatially uniform steady states i.e. fjiu^u^) = 
0, j = 1,2. To examine the linear stability of (u\,u* 2 ) let 

Uj = u* + Ej exp(ikx + At), j — 1,2. (5) 

One gets 



A + Dik 2 - an -ai2 + D 12 k 2 
-a 21 + D 21 k 2 A + D 2 k 2 - a 22 



0, 
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where 



A 



7, (°n + b 22 ) ± \/ (fen + b 22 ) 2 - 4(611622 - b 12 b 21 ) , 




9/i(«*,«2) 



Turing instability occurs when (u*, i^) is asymptotically stable if diffusion 
does not exist (k = 0); but unstable if diffusion exists (k 7^ 0). Thus the 
conditions for Turing instability are 



The third condition implies 

(D 1 D 2 - D 12 D 21 )k 4 - (D ia22 + D 2 a n - D 12 a 21 - D 21 a 12 )k 2 + 
(aua 22 — a i2 a 2 i) < 0, Vfc. 

Minimizing the left hand side of the above inequality with respect to k one 
gets the conditions for Turing instability in the following form 



{Dia 22 +D 2 au-D 12 a 21 -D 21 a 12 ) > 2 J (a u a 22 - ai 2 a 2 i)(D 1 D 2 - D 12 D 21 ) > 0. 



Now we consider Turing instability for the system (2). Following the 
above procedure one gets: 

Proposition 1: The system (2) is Turing stable if 



Proof: For the system (2), we have 

an = -fi\a 2 < 0, a 12 = -a 2 < 0, a 21 = a x - (3ia 2 > 0, a 22 = 0. 
Since D 1 >0,D 2 > 0, then 



an + ^22 < 0, ana 22 — ai 2 a 2 i > 0, bub 22 — 



b\ 2 b 2 \ < 0. 




(6) 



D 21 > 0, D 12 > 



(7) 



(Dia 22 + -D 2 an - -Di 2 a 2 i - D 21 a 12 ) < 0, 
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hence the condition (6) is not satisfied. 

Now the persistence of the system (2) is discussed: 

Definition 1: A partial differential equation 

F(t, x, u(t, x), du/dt, du/dx, ...) = 0, 

defined on a spatial domain Q is persistent in the context of zero Dirichlet 
boundary conditions if there exist a region V = {u : < u < w} in Q such 
that all solutions with nontrivial, nonnegative initial conditions are attracted 
to V. 

Equivalently the boundary of the positive cone in the space where the solu- 
tions exist act as a repeller. 

Persistence [Hofbauer and Sigmund 1998] is an important concept in pop- 
ulation biology since it means that a certain species will continue to exist. 
The work of Cantrell et al [Cantrell et al 1993] is basic in finding sufficient 
conditions for persistence in some partial differential equations without cross 
diffusion. Therefore we will use it to derive the necessary conditions for 
persistence of the system (2). Here zero Dirichlit boundary conditions are 
assumed: 

Ui = at x = 0, x = L. 

For the sake of completeness we mention the following results [Cantrell 
et al 1993] with slight modifications: 

Theorem (1): 

(i) Suppose that f(~r*,u) is Lipschitz in ~r* and continuously differentiable 
in u with 

df/du < for u > 0, /( r>, u) < if u > I, 

for some constant /, and /("f > , 0) > at some point in the domain Q , 
then the problem 

du/dt = DV 2 u +uf(T,u) inttx (O.oo), (8) 
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with Dirichlet or Neumann boundary conditions has a unique positive 
steady state u ss which is a global attractor for nontrivial non-negative 
solutions (hence the system (8) is persistent) if the following problem 
has a positive eigenvalue a 

au = DV 2 u + uf(~t, 0) inQ , 

with the same boundary conditions as (8). 

(ii) Suppose that fi(1^, u±, u 2 ) and f 2 (~r*, ui, u 2 ) are C 2 ,bounded, that /i("f > , ui, 0) 
satisfies the conditions of part (i) with positive a, that f 2 (~^, 0, U2) < 
for u 2 > 0. Let er be the largest eigenvalue of the system 

a ui = DiV 2 Ui + Ui/i("f > , 0, 0) m il, 

with Dirichlet or Neumann boundary conditions. Also let o^be the 
largest eigenvalue of the system 

a 2 u 2 = D 2 V 2 u 2 + u 2 f2('r ¥ ,u s 1 s ,0) inQ, 

where u{ s is given by fi(~f^, u\ s , 0) = 0, then the following system is 
persistent if both a and a 2 are positive 

dujdt = D^ 2 ^ +«i/i(T > ,Ui,u 2 ), (9) 
du 2 /dt = D 2 V 2 u 2 + u 2 f 2 (l^, ui, u 2 ) inttx (O.oo). 

Now we consider the persistence of system (2). The first persistence 
criteria is that the prey model alone (without predators) should be persistent, 
hence [Cantrell et al 1993] 

^2 



ai>Dl \L) • (10) 
The second criteria is that the predator system should be persistent when 
the prey is at its maximum capacity U\ = a±/ /3i, hence 

!->*(i)'- 

Now we study persistence in the presence of cross diffusion. The most 
critical region is where both u\ and u 2 are small, hence we assume that 
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Uj = ej exp(At) sin ^— J , j = 1, 2, 

and linearize in Ej, we get the following sufficient conditions for persistence 
(A>0). 

Proposition 2: In addition to Eqs. (9) and (f 0) one of the following con- 
ditions is necessary for the persistence of the system (2): 



or 




Proof: Substituting Uj in (2) and linearizing in Ej one gets 

e 1 [\-a 1 + D 1 (ii/L) 2 ]+ s 2 D 12 (ir/L) 2 = 
e 2 [X + a 2 + D 2 {n/L) 2 } - £l D 21 (n/L) 2 = 

To get a nonzero solution the determinant of the coefficients has to be zero 
which gives the following equation 

A 2 + A[-(«i - a 2 ) + (£>i + L> 2 ) (f ) 2 ] + [(-«! + A (f ) 2 )(a 2 + D 2 (f )') + 
^12^21 (f) 4 = 0. 

Persistence requires that the real part of one of the roots, at least, should be 
positive hence 




This completes the proof. 
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3 Cross diffusion in epidemics 



In this section, we will study an epidemic model where the population is 
classified into susceptible (S), and infected and infective (I), according to 
the health of each individual. The transition between the two states occurs 
according to the following rules: 

Susceptible individual having at least one infected nearest neighbor is 
infected in the next time step. 

Infected individuals are recovered and become susceptible with a rate 
A. 

Cross diffusion will be included, where S-individuals diffuse away from 
I-individuals. This system will be studied both without and with delay. In 
the first case we have 

ds = _ SI 

I = S/-A/-D0, (14) 

where both A and D are nonnegative constants. 

Here we consider that only susceptibles are capable of diffusing. This 
corresponds to the case where the disease has symptoms hence susceptibles 
try to avoid infectives. 

Studying the wave solution of Eq. (14) 

S = S(z), I = I(z), z = x — ct, 
with the boundary conditions 

S(oo) ~ 1, < S(-oo) < 1, J(oo) = I(-oo) ~ 0, (15) 

one gets 

/ = -cV = -DS" -XI + SI, S' = ^L, I' = —. 

S az az 

Substituting from the first equation and integrating the second one, we get 

SI 

F-cI = -D XclnS + cS. 

c 
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To determine the constant (F), consider z — > oo i.e. the past then Eq. (15) 
implies that F = c. Now consider the future z — > — oo, and let S(t = oo) = a, 
then the fraction of survivors a is given by the transcendental equation as 
follows: 

m = x ' <16) 

which is identical to an equation obtained by Murray [Murray 1993] in an- 
other context. 

It is straightforward to see that the homogeneous solution (S — 1, 1 = 0) 
is asymptotically stable if 

A < 1, c>VD. (17) 
Introducing delay, the system (14) becomes 



9S _ or 
dt - 

91 - S(x, t - T) I(x, t — T) — XI — D® s 



(18) 

This model is similar but not identical to that of Mendez [Mendez 1998]. 
Consider the wave solution of Eq. (18) with the boundary conditions (15), 
and studying the stability of the homogeneous solution (S — 1, 1 = 0), then 
one gets 

Substituting with / = e mz , then 

(c-^jm-\ + exp(dTm)=0. (19) 

Denoting the left hand side of Eq. (19) by f(m), then it is easy to see that 
f'(m) > V c > \f~D i.e. f(m) is monotonically increasing. Furthermore 
/(0) >0V0<A<1. Thus we have 

Proposition 3: The solution (S = 1, 1 = 0) for Eq. (18) subject to 
the boundary conditions (15) is asymptotically stable for (z — > oo), if the 
conditions (17) are satisfied. 
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The derivation of the fraction of survivals a after the epidemic wave is similar 
to the case without delay and is given by Eq. (16). 

Also, we will discuss one of the mechanisms of cross diffusion namely re- 
cruitment through signals [Saffre and Deneubourg 2002]. In this case swarm- 
ing individuals e.g. insects recruit others via signals e.g. pheromone for ant 
colonies. Here we generalize the model of Saffre and Deneubourg to include 
diffusion, and study the boundedness and persistence of the solutions. 

Let u be the population density and v be the pheromone density, then 

§ = (a + ptu)u(l - u )-0uv + D^, 

ft = e uv -v + D 2 ^, (20) 

where the parameters a, [3, 6, \x are nonnegative. There are three homoge- 
neous equilibrium solutions 

<"•"> = (0 ' 0) - (1 ' 0) ' (?• ^-dm + p) )- 

The solution (0, 0) is unstable, the solution (1, 0) is asymptotically stable if 
6 < 1. The third (positive) solution is asymptotically stable if 

9 > 1, //(l -2.0) < a. 

The results of Cantrell et al [Cantrell et al 1993] are directly applicable to 
the system (20), hence the solutions of Eq. (20) are persistent if 

a>D 1 (£) 2 , 6>1 + D 2 (£) 2 . (21) 

To study the boundedness of the solutions to Eq. (20), we will consider 
only solutions satisfying the following condition 

u(0) > 0, v(0) > =>• u(t) > 0, v(t) > Vt > 0. 

It is direct to see that du/dt < V u > 1, hence sup u(t) < 1. Consequently 
dv/dt <(9-l)v =>• V 9 < 1, then v(t) < v(0). Thus we have: 
Proposition 4: 

(i) For < 9 < 1, the solutions of the system (20) with Di = D 2 = are 
bounded in R 2 + = {(x,y) such that x > 0, y > 0}. 
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(ii) The solutions of Eq. (20) are persistent given that the conditions (21) 
are valid. 



4 Cross diffusion and evolutionarily stable strat- 
egy 

Evolutionary stable strategy ESS [Hofbauer and Sigmund 1998] is an im- 
portant concept in population dynamics. Consider a population in which 
each individual adopts one of k possible strategies. A strategy I* is ESS if 
for all I ^ I* (yet close to it), then 

(i) A(I*,r) > A(I,I*); or 

(ii) If A(I*, I*) = A(I, /*), then A(I*, I) > A(I, I), 

where A is the payoff matrix. The dynamics of ESS is given by replicator 
dynamics 

^ = Vl [(Ap)i - P Ap] , i = 1, 2, 3, k, (22) 

where Pi is the fraction of the population adopting the strategy i. 

Recently Hofbauer has introduced the concept of spatial game [Hofbauer 
1999], where the system is given by 

fa= Pi [(Ap) i -pAp] + D i ^. (23) 

This concept is generalized to asymmetric games [Ahmed et al. 2002]. Here 
we generalize it further to contain cross diffusion in addition to ordinary 
diffusion. The equations become 

f ^Pri,-^] + A^ + EA,|k (24) 

3^ 

The condition J2iPi — 1 implies A — D k + J2*j^i Dij = 0. 

Applying to the hawk(H)-dove(D)- retaliate(R) game [Maynard Smith 
1982] where players either fight for a given prize with value v or agree to 
share it. The first strategy is called hawk (H) and its drawback is that fight 
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costs value c> v. The second strategy is called dove (D) and its drawback is 
that it loses against hawks. The third strategy is called retaliate where the 
player imitates the strategy of the opponent. Hence he (or she) is a hawk 
(dove) if the opponent is a hawk (dove). Thus the payoff matrix A is 



v 



y v — c 

f 
[v-c) I 



v — c) 



where c and v are constants and c > v > 0. Let Pi(p 2 ) be the fraction 
of population adopting H(D) strategy, respectively, then the system (24) 
becomes 



dpi 
dt 



dp 2 
dt 



Pi 



D 1 



= P2 



C 

2 



(Ppi 
dx 2 



Pi , P2(V + C) p\c 



+ 



D 



12 



Pi c - o 



d 2 p 2 
dx 2 ' 
p\c 



cpip 2 



dx 2 



dx 2 



(25) 



Assuming 1 > x > 0, then linear stability analysis shows that p\ = p 2 = is 
unstable if 



^12^21 > D 2 (D 1 + 



2tt 2 



(26) 



Thus cross diffusion can destabilize equilibrium solution, in some cases, if 
D u , D 2 \ are large enough. 



5 Cross diffusion and the dynamics of learn- 
ing in multiagent systems 

Recently [Sato and Crutchfield 2002] have studied the dynamics of learning 
in multiagent systems, where the agents use reinforcement learning. They 
showed that, although the agents are not directly interacting with each other, 
a collective game between them arises through their interaction with the en- 
vironment. Such interactions can be modelled via replicator type equations. 

Consider two agents u and v with n possible actions, let Uj(t) {vj(t)) be 
a measure of the probability that the agent u{y ) will take the action j at 
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time t. Sato and Crutchfield [Sato and Crutchfield 2002] have shown that 
Uj(t),Vj(t) satisfy the following equations 

^ = Uif3i [(Av)i - uAv] + OL X Ui uj In , 

civ ■ / V \ 

= v l (3 2 [(Bu) t -vBu] + a 2 v l J2vj^[-^) , (27) 

where A and B are payoff (reward) matrices for u and v, respectively, and 
«i, a 2 , Pi, are nonnegative constants. If 0:1 = 0:2 = 0, then one regains 
the standard asymmetric replicator equations [Hofbauer and Sigmund 1998, 
Ahmed et al. 2002]. For simplicity we consider the following payoff matrices 
A = diag(ai,a 2 ), B = diag(&i,& 2 ), > 0, bi > 0, hence the system (27) 
becomes (ui = u, u 2 = 1 — u, v± — v, v 2 — 1 — v) 

du 

— = - v) (u(ai + a 2 ) - a 2 ) , 

^ = + 63) -63). (28) 

or 

There are four equilibria (u, u): £0 = (0,0), i?i = (a 2 /(a 1 + a 2 ), 0) , i? 2 = 
(0, b 2 /(b 1 + b 2 )) , _E 3 = (1, 1). The solution E is asymptotically stable while 
Ei,E 2 ,E 3 are saddle points, hence 

Corollary 1: The system (28) is not persistent. 

Including cross diffusion the system becomes 

du (9 2 t> 
— = p 1 u(l-v)(u(a 1 + a 2 )-a 2 ) + D 12 —, 

^ = p 2V (l- u )( v (b 1 + b 2 )-b 2 ) + D 21 ^, (29) 

where 1 > x > 0, and the following Dirichlet boundary conditions are as- 
sumed u(0, t) = u(l, t) = v (0, 1) = v(l,t) = 0. Assuming 

u = u eq + Wo exp(cxt) sin nx, v — v eg + v exp(crt) sin nx, 

where (u eq , v eq ) is one of the equilibrium solutions E X ,E 2 , E , one gets 
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Proposition 5: The following condition is necessary for persistence of the 
system (29): 

n n f Pip2a 2 b 2 I3il3 2 aib 2 PiP 2 ha 2 1 , OAN 

^12^21 > max ^ , — ■ -, S . 30 

{ 7r 4 7r 4 (ai + a 2 ) 7T 4 (6i + 6 2 ) J 



Proof. The condition (30) implies that Ei,E 2 ,Eq are unstable which is 
a necessary condition for the persistence of the system. 



6 Conclusions 

Conditions for Turing instability and persistence in the presence of cross 
diffusion are derived. The predator-prey model is generalized to the cross 
diffusion case. Several types of interactions are assumed. . 

Stability and persistence of the solutions of an epidemic model with cross 
diffusion is studied. 

The effect of cross diffusion to evolutionary games is studied. In some 
cases, cross diffusion destabilizes the equilibrium solution. 

Conditions for persistence and stability for the dynamics of learning in 
multiagent systems are concluded. 

Acknowledgments: 

We thank R. S. Cantrell and C. Cosner for their comments. Also we 
thank the referees for their helpful comments. 

References 

Ahmed E., Hegazi A. S. and Elgazzar A. S. (2002), On spatial asym- 
metric games. Advances in complex systems 5, 433-444. 

Cantrell R. S., Cosner C. and Hutson V. (1993), Permanence in eco- 
logical systems with spatial heterogeneity. Proceeding of the Royal. 
Society of Edinburgh 123A, 533-559. 



14 



Chattopadhyay J. and Chatterjee S. (2001), Cross diffusional effects in 
a Lotka-Volterra competitive system, Nonlin. Phen. in Complex Sys. 
4, 364-369. 

Chen L., Hsiao L. and Li Y.(2003), Strong solution to a kind of cross 
diffusion parabolic system, Comm. Math. Sci. 1, 799. 

Hofbauer J. (1999), The spatially dominant equilibrium of a game. 
Annals of Operations Research 89, 233-251. 

Hofbauer J. and Sigmund K. (1998), Evolutionary games and popula- 
tion dynamics. Cambridge university press, Cambridge. 

Le D. (2003), Cross diffusion systems on n-dimensional domains, Elec- 
tronic J. Diff. Eqns, Conference 10, 193. 

Maynard Smith J. (1982), Evolution and theory of games, Cambridge 
Univ. Press U.K. 

Mendez V. (1998), Epidemic model with infected infectious period. 
Phys. Rev. E 57, 3622-3624. 

Murray J.D.. (2003), "Mathematical biology", Springer Publishers.. 

Okubo A. (1980), Diffusion and ecological problem. Springer- Verlag, 
Berlin.. 

Saffre F. and Deneubourg J. L. (2002), Swarming strategies for coop- 
erative species. J. Theor. Biol. 214, 441-451. 

Sato Y. and Crutchfield J. P. (2002), Coupled replicator equations for 
the dynamics of learning in multiagent systems, Nlin AO/0204057. 



15 



